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! Abstract 

> 



In sterile neutrino scenarios, radiative corrections induce mass splittings propor- 



CN ' tional to the top Yukawa coupling, in contrast to the three active neutrino case 

where the induced splittings are proportional to the tau Yukawa coupling. In 

ON • view of this, we have analyzed the stability of the four-neutrino schemes favored 

ON ' 

by oscillation experiments, consisting in two pairs of nearly degenerate neutrinos 
separated by the LSND gap. Requiring compatibility with the measurements of 
the abundances of primordial elements produced in Big Bang Nucleosynthesis, 
we find that when the heaviest pair corresponds to the solar neutrinos (mainly 
an admixture of v e — v s ) the natural mass splitting is 3-5 orders of magnitude 



larger than the observed one, discrediting the scenario from a theoretical point 
of view. On the contrary, the scheme where the heaviest pair corresponds to the 
atmospheric neutrinos (mainly an admixture of — v T ) is safe from radiative 
corrections due to the small sterile component of these mass eigenstates. 
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1 Introduction 



Strong indications of neutrino oscillations have been gathered in the atmospheric and 
solar neutrino experiments |],@,[|]. These data point to two distinct mass squared 
differences: they range from 5 x 10~ 4 — 10~ 2 eV 2 for the atmospheric oscillation to 
5 x 10 -11 — 2 x 10~ 4 eV 2 for the solar one. There is further indication of neutrino 
oscillation in the, so far unconfirmed, LSND data which would require a very 
distinct mass difference: 0.3 — 1 eV 2 . In order to accommodate all mass squared 
differences, at least four light species of neutrinos are needed, one of which must be 
sterile, i.e., singlet under the Standard Model (SM) gauge group 

Among all neutrino spectra involving one sterile neutrino there are two somewhat 
preferred by the combined analysis of solar, atmospheric and LSND experiments J7|. 
These schemes are characterized by two pairs of nearly degenerate neutrinos, accounting 
for the solar and atmospheric oscillations, separated by a gap of 0(1 eV), responsible of 
the LSND anomaly (see fig.(l)). In one of these (scheme A), the two heaviest neutrinos 
correspond basically to the solar neutrinos, and in the other (scheme B), these are 
mainly the atmospheric neutrinos. The atmospheric signal is interpreted as a — v T 
oscillation while the solar deficit is explained as a v e — v s oscillation. In this paper, we 
will restrict ourselves to these scenarios. It is important to point out that when the 
solar neutrino problem is interpreted as u e — v s oscillations, only the small angle MSW 
solution (SAMSW) is allowed §j, yielding the fits 3 x lfr 6 eV 2 < Am 2 sol < 1(T 5 eV 2 
and 2 x 1(T 3 < sin 2 26 sol < 2 x 1(T 2 . 

Here we will consider a scenario with a 4 x 4 mass matrix for the three active 
neutrinos plus the sterile one. We define the effective mass term for the four species of 
neutrinos in the flavor basis as 

C = — v T M v v + h.c. (1) 

where v T = (y e , u^, u T , v s ) and A4 U is the neutrino mass matrix. This is diagonalized 
in the usual way 

M v = V*DV ] (2) 
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Figure 1: Schemes for the neutrino spectrum 
where V is a unitary 4x4 matrix relating flavor eigenstates to mass eigenstates 
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and D = diag(m 1 e' 1 ^ 1 ,m 2 e* (?i2 ,m^e % ^ 3 ,7714), with all mj > 0. It will be useful in the 
analysis of the complex case to absorb the Majorana phases in the Maki-Nakagawa- 
Sakata (MNS) matrix U 0, defined as 



U — V ■ diag(e 



-ifa/2 „-i0 2 /2 „-i<t>3/2 



1) 



(4) 



so that U 7 M. V U = diag (mi, 1712,7123,1714). 

We will consider just the SM fields plus one singlet fermion, to be regarded as 
the sterile neutrino (the extension to the Minimal Supersymmetric Standard Model 
(MSSM) is straightforward). The Majorana mass matrix for the three active neutrinos 
and the sterile one is formed by three kinds of terms. The active-active entries of the 



mass matrix come from a dimension five operator |T0 

1 



5C 



where H is the SM Higgs, Lj 



}Kij(H -Li)(H -Lj) + h.c, 



(5) 



e, /1, r) are the lepton doublets and is & 3 x 3 

2 



symmetric matricial coupling. There is also a purely Majorana mass term for the sterile 
neutrino 



5C = --v s m s v s + h.c, 
and Yukawa couplings, responsible for the active-sterile mixing 



(6) 



5C = yii> s H ■ Li + h.c. 



(7) 



After electroweak symmetry breaking, a 4 x 4 mass matrix for the four species of 
neutrinos arises _ 

where y v = (y e , y^, y T ) are the neutrino Yukawa couplings, and v = 246 GeV. Note that 
all the entries of M. v must be at most of 0(1 eV) to obtain a significant mixing between 
the sterile and the active neutrinos []. The parameters appearing in the Lagrangian 
that form the neutrino mass matrix have a different renormalization group behavior, 
as can be checked directly from their renormalization group equations (RGEs) glH 



167T 2 — 

dt 



-2,g 2 2 + 2A + 2S 



1 

K 

2 



Y+Y e + yty,) + (YlY, 



y*y 



e * e 



(9) 



(fyu 

dt 



16tt 5 



\gl + \al - s) I 3 - I (yiy, - Yt Y< 



(10) 



where 



dm s 
~dT 



8tt 2 



'IV. 



S = Tr(3Y{jYu + SY^Yd + yty, + Y+Y. 



(12) 



and t = log fi. In the formulas above, #2 and A are the SU{2) gauge coupling and the 

quartic Higgs coupling respectively, and Yu,D,e are the Yukawa matrices for up quarks, 

down quarks and charged leptons. Previous papers p!3|Jl^ , p^ , p^ . p^ , p!8| , p!9] , pO| , pT| , p2| have 

considered the modification of the mixing angles and mass splittings due to the last 

term in eq.(|D|) in scenarios with three active neutrinos. The first term in eq.(^), which 

^^This condition has usually been claimed to be rather unnatural. However, see JTT[ for a framework 
derived from string theories in which this requirement can be fulfilled naturally. 
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depends on h t , produces a flavor independent modification of the active neutrino mass 
matrix. However, in eg. flT0|) the dependence on h t is different, while eq.flTTD lacks a h t 
dependent term, thus important modifications of the mass texture are expected now. 
This will have important consequences on the theoretical plausibility of some kinds 
of neutrino spectra. The aim of this paper is to impose naturalness constraints on 
sterile neutrino scenarios from radiative corrections. This approach has already been 

but 



employed for the three active neutrinos scenario ||i4|14|15| . |iq , ^ 
now the effects are expected to be larger, as a consequence of the different nature of 
the sterile neutrino. Furthermore, similar strong renormalization effects are expected 
in more general models containing sterile neutrinos. 

Recalling that the Yukawa coupling between the active and the sterile neutrinos 
must be small, in order to have masses of 0(1 eV), in a first approximation we can 
neglect all couplings but h t in eqs.(P, [IO| , |TT|) . In this approximation, the RGE for the 
neutrino mass matrix is 



dM l 
dt 



1 



16vr 2 



M v 



Sh 2 t Is 




+ 



3h 2 t l 3 











M, 



(13) 



It is worth noticing that the RGE for the mass matrix has formally the same 
structure as the one considered in previous works for the three neutrino case (see 
eq.(14) in p2| ) with P = —j^s (^r^ 3 an< ^ K u = 0- Therefore, we can apply the 



16th y o 0, 

general results derived there to obtain the RGEs for the mass eigenvalues 



drrii 
~dT 



1h 2 

B(l-I^| 2 ) m h 
8tt z 



(14) 



and the MNS matrix 



dU 
~dT 



UT, 



(15) 



with T a 4 x 4 matrix defined as 



T a = 0, 



T ■ 



167T 2 



rrii + rrij 



rrii 



rrij 



rrij — mi 



ImiU^U, 



S3 1 



^ ^ J, (16) 



rrii + m j 

with i = 1,2,3,4. Starting with degenerate neutrino masses at a high energy scale A, 
the degeneracy will be lifted at Mz by terms proportional to the top Yukawa coupling, 
in a way that depends on the last row of the MNS matrix. Strictly speaking, when 
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rrii = rrij, there would be a divergence in eq.(jI3[), but there is also an ambiguity in 
the definition of the MNS matrix. It can be shown |22j that, once the ambiguity is 
removed by the running, the mixing matrix satisfies RelU^Usj) = 0, making eq. (|T5|) 
finite. 



For the MSSM case, the relevant RGEs for this analysis can be found in ||16| . |23 
although there is a complete analogy with the SM case: neglecting all couplings but 
the top Yukawa coupling, eq. flliD is still valid with the replacement h t — > h t / sinfi. 

In view of the form of the RGEs for the masses and mixing angles, any information 
about the U S i elements of the MNS matrix would be certainly valuable. Fortunately, 
the remarkable agreement between the predictions of Big Bang Nucleosynthesis (BBN) 
and the observed abundance of primordial light elements sets stringent bounds on U S i. 
The upper bound N BBN < 4 for the effective number of neutrinos implies that, for 
scheme A 

|^ 3 | 2 + |^4| 2 ^0.99, (17) 

while for scheme B 

|f/ s3 | 2 + |^ 4 | 2 ^10- 4 . (18) 
These constraints will have important consequences in the analysis to follow. 



2 The effect of radiative corrections 

As explained in the introduction, radiative corrections play an important role when 
two mass eigenstates are almost degenerate, i.e., Am| <C m?j. If this is the case, the 
renormalization group evolution can produce mass differences larger than the observed 
ones, requiring a fine tuning of the neutrino parameters at high energy to achieve a 
correct splitting. This fact was exploited in pTT|JT8| ] to impose severe constraints on 
the vacuum oscillation solution to the solar neutrino problem. In this paper we are 
going to perform a similar analysis for the experimentally favored neutrino patterns 
presented in the introduction. 

We will assume that at some high energy scale A new physics generates a neutrino 
mass texture A4 U (A) in which the heaviest pair of neutrinos have masses nearly degen- 
erated and the other two are ~ (we neglect the effects of the RGEs on the masses 
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of the two lightest eigenstates because a hierarchical spectrum will not be essentially 
modified by the running). After integration from A to Mz, the radiatively corrected 
mass eigenvalues can be obtained from eq. flT3|) , giving 

m 1 {M z ) ~ 0, 
m 2 (M z )~ 0, 

(19) 

m 3 (M z ) ~ m 3 (A) J t J s3 , 
m 4 (M z ) ~ 7714(A) 7 t J s4 , 

where 

The typical size of I t is ~ 0.93 (~ 0.67) for A = 10 3 GeV (10 12 GeV) in the SM, and 
~ 0.52 for A = 10 12 GeV in the MSSM with tan/? = 2. 

In the following subsections, we will analyze two different patterns, labeled as 
schemes A and B. In scheme A, the heaviest neutrinos participate in solar oscilla- 
tions, and thus the mass splitting between them is small. In addition, these mass 
eigenstates have a large sterile component. Therefore, it is in this scheme where larger 
renormalization effects are expected and will be studied in detail, both the real and 
the complex case. On the other hand, in scheme B the heaviest pair of neutrinos has 
a small sterile component and the required splitting between them is larger, so in this 
case radiative corrections are not expected to have dramatic consequences. 

2.1 Scheme A 
2.1.1 Real Case 

Throughout the analysis of the real case we will work for convenience with the V matrix, 
setting all the CP violating phases in V to and allowing the Majorana phases to be 
either or n, yielding positive or negative mass eigenvalues, respectively. In this case 
the renormalization effects are sensitive to the relative sign of the mass eigenvalues. In 
consequence, we will consider the two possibilities separately: 7713 ~ 7774 and 777.3 ~ —7774. 

a) 7773 ~ 7774 

The solar splitting produced in scheme A by the renormalization group evolution 
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can be obtained from eqs.([L9|), giving 



AmU-< ^(I-^Vm), (21) 

where m\l 2 l 2 ^ = m\(Mz) — Am 2 LSND ~ 1 eV to account for the LSND anomaly. 

The degeneracy of the mass eigenstates at A has also an important impact on the 
radiative corrections of the mixing matrix. The renormalization group equation of the 
MNS matrix, given by eq . (|l5l) , is dominated in this case by 



lo7r z m 4 — 777,3 



Thus, the renormalization group evolution will quickly drive V S 3 — > in the infrared 
22|. Therefore, the BBN constraint fllTD implies that V£(M Z ) £ 0.99. This value is 



close to the fixed point V S i(Mz) = 1 and will not change appreciably with the scale, 
hence, it can be considered constant. In this approximation, and neglecting terms of 
the order of 1 — V 2 4 , eg. (|2"ID reads 

AmL * Am 2 LSND (l - f t Vh ) ~ AmU(l - (23) 

For example, taking Am 2 LSND ~ 0.3 eV 2 the induced splitting is Am^, ~ 0.04 eV 2 (0.17 
eV 2 ) for A ~ 10 3 GeV (10 12 GeV) in the SM, and Am 2 sol ~ 0.22 eV 2 for A ~ 10 12 GeV 
in the MSSM with tan (3 = 2, which are too large compared to any explanation of the 
solar neutrino deficit by several orders of magnitude. 

The only way out of this shortcoming would be an extreme fine tuning between the 
initial values of the mass splittings and the effect of the RGEs. To obtain the correct 
value for Am 2 ol at M%, the mass splitting at the scale A is forced to lie in the narrow 
range 

|m 2 (A) - m 2 (A)| ~ \{I; 2 - l)Am 2 LSND ± Am 2 J, (24) 

where we have neglected again terms of the order of 1 — V, 2 . In the most favorable 
choice of scenario and parameters for the fine-tuning problem, namely the SM with 
A = 10 3 GeV and Am 2 LSND = 0.3eV 2 , we get from (||) |m 2 (A) - m 2 (A)| ~ (2.309 ± 
0.001) x 10 -2 eV 2 , i.e. the initial splitting has to be adjusted in less than one part 
in 10 3 . Raising the values of Am 2 LSND or A, or going to the supersymmetric case, 
requires a higher degree of fine tuning. This subtle conspiracy between initial masses 
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and radiative corrections makes this scenario implausible from the point of view of 
naturalness. 



b) m 3 ~ — m 4 

In this case, the neutrino mass eigenvalues at M z are given by eq.([L9|), with m 3 (A) = 
—777,4 (A) = m o, while the solar mass splitting induced by the radiative corrections is 
given by eq.(PTl). Now we do not expect abrupt changes in the relevant elements of the 



mixing matrix for the solar splitting (there is no ambiguity in the diagonalization of 
the mass matrix). To be more precise, the RGEs for these matrix elements read 

iV ' 3 ~ :ik ':V, 3 (l-Vi-V^ (25) 



dt 167T 2 

(for V S 4 the RGE is similar with the interchange 3 <-> 4). The BBN bound (|T7| ) implies 



-v 2 

that these elements can be considered constant in first approximation, hence I S i — It 



It is apparent from (pi]) that the splitting will be unacceptably large unless V S 3 ~ V^. 

For definiteness, we will use the simple texture for the mixing matrix given in 
|[25| , although our results will not change essentially if we allow for more complicated 
textures. This texture almost decouples the v^v s — v\, u 2 and the v T — u 3 , z/ 4 sectors 
(as suggested by BBN bounds) and slightly rotates the v\ — v% plane by a small angle 
V^Olsnd to generate a small mixing between the v e ^v s — ui,u 2 and v T — v 3 , z/ 4 
sectors, in order to accommodate the LSND signal. We could also have small rotations 
in the v\ — 7/4, v 2 — V3 and z/2 — planes, but, as stressed before, the main results will 
remain unchanged. The mixing matrix at Mz is 



V{M Z 



/ cos 9 sol sin V29 LSND cos 9 so i cos V28 LSND sm9 sol \ 

COS 9 atm COS \f29 LS ND sm9 at m - cos 9 atm sin y/29 L SND 

- sin 9 at m cos V29 LS nd cos9 atm sin 9 atm sin V29 L snd 

V sin 9 so i sin V29 LS nd sin 9 so i cos V^9 LS nd cos 9 soi ) 



(26) 



Using this mixing matrix and eq. (|21|) , it is easy to see that the renormalization group 
evolution produces a solar mass splitting given by 

A™ 2 ~ A™ 2 M T -2(cos 2 8 sol -sin 2 8 sol cos 2 V2e LS ND)] ( 07 \ 

Lym sol — ^ m LSND[ L ~ 1 t \ \ zt ) 

It is obvious that in order to obtain a sufficiently small splitting we must require 
cos 2 9 so i ~ sin 2 9 soi cos 2 V^9 LSND (typically, to obtain a mass splitting of the solar order 
the condition | cos 2 9 so i — sin 2 9 so i cos 2 V%9lsnd\ < must be fulfilled). Due to the 
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smallness of Olsnd we infer that the mixing between v e and u s must be nearly maximal, 
and adjusted with high accuracy. This discredits the small angle MSW solution to the 
solar neutrino problem (the only acceptable one) from the point of view of naturalness. 
The only way out would be, as before, an extremely artificial fine tuning between the 
initial masses and the renormalization group effects. 

2.1.2 Complex Case 

In the introduction it was explained that when m 4 (A) = m 3 (A) the running quickly 
drives the mixing matrix to a form in which Re(U* 3 U S 4) = 0. In this section we will 
assume that we end up with a phenomenologically viable mixing matrix and then we 
will investigate the effect of the running on the mass splitting. 

When Re(U* 3 U S 4) = 0, the RGEs for \U s3 \ 2 and \U si \ 2 can be written as 

-^sl 2 - -ff (1 - \U s3 \ 2 - |^4| 2 )|t/ S 3| 2 + ^|^^|^3| 2 |f/ s4 | 2 , (28) 
at 87T Z 6ir z m 4 + m 3 

(for \U S 4\ 2 the RGE is the same with the interchange 3 «-> 4). After imposing the BBN 
constraint ([17|) we can approximate the above equation by 

Lvj * -'\vj = 3 A"i luMUsi f. (29) 

at at Sir z 1714 + m 3 

As said before, in the four-neutrino scenarios considered here only the small angle MSW 
solution is allowed |§, then \U s3 \ 2 is small and by the BBN constraint (|i~7|) , |£/ S 4| 2 is 
close to one. Hence, eq.(p9|) implies that both \U s3 \ 2 and \U s4 \ 2 are approximately 
constant with the scale, yielding I s3 ~ I t ' (7s3 ' and J S 4 ~ J"'* 754 ' . Under those well 
founded hypotheses, the solar splitting induced by the running is 

AmL - Am 2 LSND (l - I 2 ^-^) „ Am 2 LSND (l - I 2 ), (30) 

where we have neglected terms of the order of \U s3 \ 2 and 1 — |?7 S 4| 2 . The numerical 
results are identical to those after eq.(^), again too large by 3-5 orders of magnitude. 
We conclude that it is impossible to reconcile a phenomenologically acceptable MNS 
matrix and an experimentally allowed mass splitting without a fine tuning between 
initial conditions and renormalization group effects. This permits us to completely 
discredit the scheme A for the four neutrino spectrum from a theoretical point of view. 
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2.2 Scheme B 



Things work quite differently in this case. For the SM, the induced splitting between 
the heaviest pair, to be identified with the atmospheric oscillation, is given by 



Am atm ~ 



r2 r-2x 
J s3 J s4 li 



(31) 



with m\I 2 I\ 



m 2 (M z 



/S.m 2 LSND . The BBN constraint (|1 



requires both \U. 



,s3 1 



and \U S 4\ 2 to be smaller than 10 4 . This condition is approximately maintained along 
the running, as can be checked from 



dt V'~ s31 ' ) - 8tt 2 

This allows the following expansion of eq.(|31 



+ \U : 



3h 2 



(1 



Am atm ~ 



U S 3\ 2 — \UgA 



\u. 



,s3 



2 + |^4| 2 ). 



2Ami SiVD (|f/ s4 | 2 



,s3| 



lOK It- 



(32) 



(33) 



Then, the atmospheric splitting induced by the RGEs is ^ 2 x 10 _4 Am| /5ArD , too low 
to account for the atmospheric signal. Inversely, a mass difference of the order of the 
atmospheric oscillation generated at a high scale A will not be upset by the radiative 
corrections. 

For the MSSM case, due to the smallness of |?7 S 3| 2 and |t/ S 4| 2 , the effects induced 
by the top Yukawa coupling will also be small. However, the effects of the tau Yukawa 



coupling can be dominant in the large tan/3 regime ||15| , |16| , so it is worth considering 
them in more detail. After including this coupling, the renormalization group equations 
for the mass eigenvalues are 



drrii 



3h* t 



:i-\u s , 



+ 



hi 



\U Ti \ 



rrii 



(34) 



_87r 2 sin 2 /^ '~ sl1 ' ' 87r 2 cos 2 /3' 

where i = 1,2,3,4 and h T is the tau Yukawa coupling. The induced splitting, in the 
leading-log approximation, is given by 



^■ m atm — 



2Am 2 LSND 



k t {\U : 



si 



+ k T (\U rZ 



\Ur4\ 2 ) 



log(A/M 2 



where 



3h 2 



k T 



(35) 



(36) 



87r 2 sin 2 /5' '" T 8tt 2 cos 2 /5' 

2 , \U S 4\ 2 and h T are smal 
In fig. (2) it is plotted the maximum splitting reachable by the RG evolution for A = 10 12 



(The leading-log approximation is well justified because |t/ S 3| 2 , | C^ S 4 1 2 and h T are small.) 
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4xl0" 3 




Figure 2: Atmospheric splitting for A = 10 12 GcV and A = 10 16 GeV. The dashed line corresponds 
to the lower bound for Am 2 tm . 

GeV and A = 10 16 GeV. We have assumed the most favorable case, i.e., Am 2 LSND = 
1 eV 2 , (\U s4 \ 2 - \U s3 \ 2 ) ~ 10~ 4 and sin 2 6 atm ~ 0.82, so that (\U r3 \ 2 - \U r4 \ 2 ) ~ 0.42. 
From fig. (2) it is apparent that in order to generate the atmospheric splitting radia- 
tively, a moderately large tan/3 is needed, or a sizeable value of the cut-off. 

3 Conclusions 

We have shown that radiative corrections in models with sterile neutrinos are poten- 
tially dangerous in many scenarios. Due to the different interactions of the sterile 
neutrino, the renormalization group behavior of the mass matrix for the active and 
sterile neutrinos presents a strong non-universality. This leads to induced splittings 
proportional to the top Yukawa coupling, in contrast with the three neutrino scenario, 
where these splittings were found to be proportional to the tau Yukawa coupling. This 
strong scale dependence of the mass matrix texture has allowed us to make some 
statements concerning the viability of some neutrino spectra. More precisely, we have 
analyzed the two schemes favored by the oscillation data, both consisting in two pairs 
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of nearly degenerate neutrinos separated by the LSND gap. Requiring compatibility 
with the measurements of the abundances of primordial elements produced in Big Bang 
Nucleosynthesis, we have found that when the heaviest pair corresponds to the solar 
neutrinos (mainly an admixture of v e — u s ) the splitting is not stable under radiative 
corrections, and the induced one is 3-5 orders of magnitude larger than the observed 
splitting, making the scenario completely implausible from a theoretical point of view, 
the only way out being an extreme fine tuning in the high energy parameters. Inciden- 
tally, this scenario can be tested in the near future with neutrinoless double-/? decay 
experiments ||26|| . 

On the other hand, when the heaviest pair corresponds to the atmospheric neutrinos 
(mainly — u T ), the scenario is safe from radiative corrections due to the small mixing 
of these mass eigenstates with the sterile neutrino. For the SM, the induced mass 
differences are lower than the atmospheric splitting (although close to the lower bound 
for this oscillation) while for the MSSM, there are regions of the parameter space 
where, beginning with degenerate neutrinos, the atmospheric splitting can be generated 
radiatively for a moderately large tan (3 or a sizeable cut-off. 
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